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The invariant Buda-Lund parameterization of the two-particle Bose-Einstein cor- 
relation functions is presented, its derivation is summarized. In its particular 
multi-variate Gaussian limiting case, the invariant Buda - Lund parameterization is 
compared to the Bertsch-Pratt and Yano-Koonin-Podgoretskii parameterizations, 
advantages and shortcomings are discussed. The invariant Buda-Lund parame- 
terizations are given also for non-Gaussian multi-variate distributions, including 
damped oscillations in the like-particle correlation function, that are similar to the 
oscillating intensity correlations of binary stars in stellar interferometry. A separa- 
tion between the pion and the proton source is also estimated in the Buda - Lund 
hydro framework, the result is utilized to extract the mean proper time of particle 
emission with the help of fits to E877 data on non-identical particle correlations 
by Miskowiec. 



I Introduction 

Presently there is an increasing interest both in high energy heavy ion and in 
high energy particle physics to describe in greater details the space-time picture 
of particle emission with the help of particle interferometry: the space-time 
geometry of the freeze-out hypersurface carries information on the existence of 
a transient Quark-Gluon Plasma phase to be created in heavy ion collisions. 
Overlapping regions of (qq) + (qq) jets from a decaying W + , W~ pair at LEP- 

II may result in unwanted systematic errors in precision determination of W 
mass in particle physics, and the magnitude of this effect can be estimated 
with precision only if the space-time picture of these reactions is reconstructed 
and the magnitude of Bose-Einstein correlations between pions from different 
W-s is calculated correspondingly. 

An invariant formulation of Bose-EinsteirL£orrelation functions was found 
by the Budapest-Lund collaboration in refs. EH3. This parameterization is re- 
ferred as the Buda-Lund parameterization of the correlation function, or BL in 
short. Although the BL results are rather generic, they contain in particular 
limiting cases the power-law, the exponential, the double-Gaussian, the Gaus- 
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sian or other parameterizations. The BL parameterization of the correlation 
function yields not only a boost-invariant description, but also a rather simple 
functional form. 

We consider high energy heavy ion reactions or single jets in high energy 
particle physics coming from a fragmentation of an energetic leading quark. 
In these physical situations, a dominant direction of expansion of the parti- 
cle source is identified. We denote space-time coordinates by x — (i, r) = 
(t,r Xl r y ,r z ) and momentum variables as k = (E, k) = (E,k Xl k yi k z ). We 
choose the direction labelled with subscript z to coincide with the dominant 
direction of the expansion. 

The single-particle spectra and the two-particle correlations are deter- 
mined in the Wigner-function formalism. As this formalism is fairly well-known 
now, we move the summary of the derivation to section ||, and jump to the 
results immediately, as given in the next section. We then compare the BL 
form of the correlation function to the Bertsch - Pratt (BP) and the Yano- 
Koonin - Podgoretskii (YKP) 3-dimensional parameterizations in a particular, 
Gaussian limiting case of the BL correlation function, as the BP and the YKP 
parameterizations are defined at the moment only for Gaussian correlation 
functions, as far as we know. We also considerj-effective separation of sources 
in the BL hydro model, also described in ref. El. The effective separation of 
like-particle sources is shown to result in small, damped oscillations in the in- 
variant longitudinal component of the two-pion intensity correlation function. 
The effective separation of unlike-particle sources in the BL hydro parameteri- 
zation is larger. Such a separation was recently determined experimentally by 
D. Miskowiec et al in the analysis of E877 data on Au + Au collisions at BNL 
AGS, based on their estimated L = 10 fm/c separation scale and the BL hydro 
picture we estimate the mean freeze-out time in Au + Au collisions. 



2 The Buda-Lund parameterization for Bose-Einstein correlations 

The two-particle correlation function is defined as 

, iV 2 (ki,k 2 ) , . 

C(ki,k 2 ) - — — — — — — , 11) 

Ai(ki)A 2 (k 2 ) 

a ratio of the two-particle invariant momentum distribution to the product of 
the single-particle invariant momentum distributions. Hence the correlation 
function is invariant by definition, see ref. □ for an exploration of its properties 
based on its Lorentz invariance. We assume, that with the help of certain ex- 
perimental methods all non-Bose-Einstein related correlations can be removed 
from this function, and we focus on the Bose-Einstein correlations only, as 
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these correlations carry information about the space-time distribution of par- 
ticle emission. The mean and the relative four-momentum are introduced as 

K = (k 1 + k 2 )/2, (2) 
Ak = k 1 -k 2 . (3) 

In the following, the relative momentum four-vector shall be denoted also as 
Afc = Q = {Qo,Qx,Q y ,Qz) = (QojQ); and in general a shall denote the 
vector part of a four- vector a — (ciq, a). 

A simple invariant formulation of Bose-Einstein correlation functions was 
given, for cylindrically symmetric, longitudinally expanding particle sources in 
refs.lM Such sources may be characterized not only with longitudinal but also 
by transversal and temporal inhomogeneities. In a Gaussian approximation, 
the Buda-Lund form of the Bose-Einstein correlation function reads as follows: 

C( Ak, K) =1+A cxp (-R 2 = Q 2 = -R\Q\~ R±Q±) , (4) 

where the fit parameter A measures a strength of the correlation function. 
The fit parameter R = reads as i?-timelike, and this variable measures a width 
of the proper-time distribution. The fit parameter i?y reads as i?-parallel, it 
measures an invariant length parallel to the direction of the expansion. The fit 
parameter R± reads as i?-perpedicular or i?-perp. For cylindrically symmetric 
sources, R± measures a transversal rms radius of the particle emitting source. 

The invariant time-like, longitudinal and transverse relative momenta are 
defined with the help of another, to this point suppressed fit parameter, rj, 
which characterizes the direction of the center of particle emission, x — (t, f z ) in 
the (t, r z ) longitudinal coordinate space. Such a direction can be characterized 

by a normalized four - vector n(x) — x/t, where n ■ n = + 1, and r = \Jt 2 — r 2 z 
is the mean longitudinal proper-time of particle emission. The such a direction- 
pointing normal vector n can be parameterized as n — (coshfTj], 0, 0, sinh^]), 
where r] = 0.5 log [(t + r z )/(t — r z )] is the space-time rapidity of the center of 
particle emission, see Fig. 1. Space-time rapidity is a space-time coordinate, 
that transforms additively in case of longitudinal Lorentz-boosts, similary to 
the the rapidity variable in momentum space. A boost - invariant decomposi- 
tion of the relative momenta can be defined, as follows. The invariant time-like, 



parallel and perp relative momentum components read as 

Q= = Qo cosh[rj] - Q z smh[rj] = Q -n (5) 

Q\\ = Qo sinh[r/] - Q z cosh[^] = Q x n, (6) 

Q± = \f0l + Q 2 y , (7) 

Q 2 = Q-Q = (Q -n) 2 - (Q xn) 2 - Q\. (8) 
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In the above, a-b = a^b^ — aobo — ab = aobo — a x b x — a y b y — a z b z stands for the 
inner product of four- vectors. As fi ■ n = +1, this direction is time- like, hence 
Q ■ n — Q= is an invariant time-like component of the relative momentum. 
Q\\ is an invariant relative momentum component in the longitudinal direction 
(parallel to the beam or to the thrust axis), and Q± is the remaining perp or 
transverse component of the relative momentum, also invariant for longitudinal 
boosts. 

In short, we have the following 5 free parameters for cylindrically sym- 
metric, longitudinally expanding sources: A, R =1 R\\, R± and rj. In principle, 
each of these parameters may depend on the mean momentum K. At any 
fixed value of the mean momentum K, the 5 free parameters of the invari- 
ant Buda-Lund correlation function can be fitted to data; alternatively they 
can be theoretically evaluated from an assumed shape of the emission function 
S(x, k), for example see rem. In the Buda-Lund parameterization, the explicit 
mean momentum dependence of the parameters can be written as follows: 

C(Ak,K) = l + A(K)exp(- J Ri(K)Qi(K)-E|(K)Q|(K)-fli(K)Qi) , 

(9) 

Note that the mean momentum dependence of the relative momentum com- 
ponents Q = (K) = Q = (n(K)) and (K) = Q||(n(K)) is induced by the mean 
momentum dependence of the direction pointing normal vector n(K), simi- 
larly to the dependence of the side and the out components of the relative 
momentum on the direction of the mean momentum in the Bertsch - Pratt pa- 
rameterization. Physically, i? = , R\\ and R± are longitudinally boost-invariant 
lengths of homogeneityca in the time-like, the longitudinal and the transverse 
directions. Hence, Rn (K) is in general less than the total longitudinal exten- 
sion of the source. Similar statements hold for the other, invariant lengths 
of homogeneity in the transverse and in the temporal directions, R±(K.) and 
R=(K). 

2. 1 Symbolic notation for the side and out components 

Up to this point, we assumed a cylindrically symmetric source, where the 
spatial information about the source distribution in (r x ,r y ) was combined to 
a single perp radius parameter R±. 

In order to distinquish easily the zero-th component of the relative mo- 
mentum Qq from the out component of the relative momentum Q D = Q ou t, 
Qo 7^ Qoj w e introduce the following symbolic notation in transverse directions: 

Qside = Q.. , (10) 
Qout = Q: , (11) 
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Q\ — Qlide + Qout — Q 2 . + Q: 



(12) 



The idea behind this notation is similar to that of symbolizing the perp direc- 
tion by index "_L" , standing for a direction that is transverse to the longitudinal 
direction. The longitudinal components of the relative momenta are symbol- 
ized by index "||", i.e. two vertical parallel lines, imagining that the z axis is 
in the plane of the paper, pointing upwards, and the considered component is 
parallel to that direction. The invariant time-like components are indexed by 
" = ", not to be confused by the equality sign =. This symbol " = " is obtained 
by a 90° rotation of the symbol "||". The remaining two orthogonal directions 
are the "side" and the "out' , we symbolize them by parallel lines as well. But 
the two possible orientations of parallel lines in the plane are used up by the 
symbols = and ||, hence these two lines are thought to be orthogonal to the 
plane of the paper, thus they are symbolized by two dots. The side component 
is described by putting the dots side-by-side, which yields index The re- 
maining out direction is orthogonal to all the above directions, we rotate side 
by 90° to obtain the symbol This symbolic notation stands for the lon- 
gitudinally boost-invariant decomposition of the relative momenta. Similarly, 
the side and the out radii can be denoted as 

Rside — R.., (13) 

Rout=R:- (14) 

Cross-terms, if any, can be denoted by straightforward mixing of the symbols, 
e.g. a possible side-out cross-term may be denoted by ". :" , an out-long cross- 
term by ": |" and a side-long cross-term by "..|" etc. In a general Gaussian 
form, suitable for studying opacity effects, the Buda-Lund invariant BECF 
can be denoted as 

C(Q, K) = l + K cxp (-RLQl - R\Q\ - R 2 .Q 2 - i? 2 Q 2 ) (15) 

Note, that this equation is identical to eq. (44) of ref. 0, rewritten into the 
new, symbolic notation of the Lorentz-invariant directional decomposition. 

The above equation may be relevant for a study of expanding shells as 
well as opacity effects as recently suggested by H. HeiselbergEI. The lack of 
transparency in the source may result in an effective source function, that 
looks like a crescent in the side-out reference frame □. The overall cylindrical 
symmetry of particle production is maintained for simultaneous rotations in the 
x and the k space, but the emission function S(x, K) at a given fixed direction 
of the mean transverse momentum K± — (0, K x , K y , 0) has different width in 
the out direction and another, possibly larger width in the side direction!!!. 
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3 General derivation of Buda-Lund shape not only for Gaussians 



First we define the correlation function with the help of the Wigner-function 
formalism, then we introduce the intercept parameter A* in the core-halo pic- 
ture. Then we evaluate the correlation function in terms of longitudinally 
boost-invariant variables and we end up with a general Buda-Lund form of 
the correlation function. 



3.1 Wigner Function Formalism 

The two-particle inclusive correlation function is defined and approximately 
expressed in the Wigner function formalism as 



c(Ak-K) = at f ;\ %- s ~ i + l : /n \J ' . (16) 



N 2 (kM „ \S(Ak,K) 
N\ (ki ) Ni (k 2 ) ~ \S(0,K) 

In the above line, the Wigner-function formalisn00'@ is utilized assuming 
fully chaotic (thcrmalizcd) particle emission. The covariant Wigner-transform 
of the source density matrix, S(x, k) is a quantum-mechanical analogue of 
the classical probability that a boson is produced at a given x^ — (t, r) = 
(t, r x ,r y ,r z ) with fc M = (E, k) = (E, k x ,k y ,k z ). The auxiliary quantity 

S(Ak,K) = J d 4 xS(x, K) exp(iAfc • x) 

appears in the definition of the BECF, with Afc = k 1 — k 2 and K = {k\ + fc 2 )/2. 
The single- and two-particle inclusive momentum distributions (IMD-s) are 
defined in terms of the cross-section, they can be evaluated as 

7v l(k ) = -^=5(Afc = 0,fc) and 7V 2 (k 1; k 2 ) = -^^-,(17) 
ot ak at aki clk.2 

where a t is the total inelastic cross-section. Note that in this work we utilize 

d 3 k 

the absolute normalization of the emission function: j — d xS(x,k) = (n). 

It has been shown recently in relk0@@0 that for explicitely symmetrized 
n-particle system of bosons with variable number of bosons, the two-particle 
Bose-Einstein correlation function is properly defined by eq. ([l6|). 

3.2 Effects from Large Halo of Long- Lived Resonances 

If the bosons originate from a core which is surrounded by a halo of long-lived 
resonances, the IMD and the BECF can be calculated in a straightforward 
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manner. The detailed description is given in refsli3'E3, here we recapitulate 
only the basic idea along the lines of ref. E3. If the emission function can 
be approximately divided into two parts, representing the core and the halo, 
S(x; K) = S c (x; K) + Sh(x; K) and if the halo is characterized by large length- 
scales so that ShiQmin] K) << S c {Qmin\ K) at a finite experimental resolution 
of Qmin > 10 MeV, then 



NiQl) =iVi, c (k)+iVi, fe (k), 

v ; |5 C (0,^)| 2 
where N± t i(k) stands for the IMD of the halo or core for i = h,c and 

^i, c (k) n2 



(18) 
(19) 



A* = \*{K = k) 



Ni(k) 



(20) 



The phcnomenological A* parameter has been introduced to the literature by 
Deutschmann long time ag<£3. In the core/halo picture, this effective intercept 
parameter A(k) can be interpreted as the momentum dependent square of the 
ratio of the IMD of the core to the IMD of all particles emitted, assuming 
completely incoherent emission from the source. 

The validity of the core/halo picture to any given reaction is not guaran- 
teed, hence systematic checks of the applicability of this simplifying picture 
has to be performed, similarly to the ones done in ref.E3. 



3.3 Invariant, generic decomposition of the Bose-Einstein correlation func- 
tion to a Buda - Lund form 

For systems expanding relativistically in one direction (r z ) , it is advantageous 
to introduce the longitudinally boost invariant variable r and the space-time 
rapidity r\ that transforms additively under longitudinal boosts, 

r = y/W^, (21) 
r] = Q.5log[{t + r z )/(t-r z )}. (22) 

Similarly, in momentum space one introduces the transverse mass m t and the 
rapidity y as 



m t = y/EZ-pl (23) 
y = 0.5\og[(E+p 2 )/(E- Pz )]. (24) 
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Fig. 1. Space-time picture particle emission for a given fixed mean 
momentum of the pair. The mean value of the proper-time and the 
space-time rapidity distributions is denoted by r and rj. As the rapidity 
of the produced particles changes from the target rapidity to the projec- 
tile rapidity the {T(y),rj(y) variables scan the surface of mean particle 
production in the (t, r z ) plane. 

In order to obtain (at least approximately) a boost-invariant picture, we 
characterize the source of particles in the boost invariant variables t, wit and 
rj — y. For systems that are only approximately boost-invariant, the emission 
function may also depend on the deviation from mid-rapidity yo. The scale 
on which the approximate boost-invariance breaks down is denoted by Arj 7 a 
parameter is related to the width of the rapidity distribution. 

A simple generalization to non-Gaussian correlators in the Buda-Lund pic- 
ture is obtained if we assume that the emission function factorizes as a product 
of an effective proper-time distribution—a space-time rapidity distribution and 
a transverse coordinate distribution Eim as 

S c (x, K)d i x = H^(r)G^(r])I^(r x ,r y ) d,TTdridr x dr y . (25) 

In the above, subscript * stands for a dependence on the mean momentum, 
the mid-rapidity and the scale of violation boost-invariance (K, y$, A77). The 
function H*(t) stands for an effective proper-time distribution (that includes 
an extra factor r from the Jacobian d 4 x = drr drj, dr x dr y , in order to simplify 
the results). The effective (K dependent) space-time rapidity distribution is 
denoted by G*(ry), while the effective transverse distribution is denoted by 
I*{r x ,r y ) . In the above equation, the mean proper-time f is factored out to 
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keep the dimensionless nature of the distribution functions. Such a pattern of 
particle production is visualized on Fig. 1. 

In ,caae of hydrodynamical models, as well as in case of a decaying Lund 
stringsoB, production of particles with a given momentum rapidity y is limited 
to a narrow region in space-time around rj and r . If the sizes of such an effective 
source are sufficiently small (or with other words if the Bose-Einstein correla- 
tion function is sufficiently broad), the plane- waves that induce the Fourier - 
transformation in the correlation function can be decomposed in the shaded 
region on Fig. 1 as follows: 

exp[«A_B(ti - 1 2 ) - iQ z (r z ,i - r z , 2 )} ^ exp[iQ = (ri -r 2 ) - iQff(r)i -Tfe)], (26) 



exp[-iQ x (r x> i - r x>2 ) - iQ v {r Vt i - r„ >2 )] = 

= ex P HQ : (r : ,i - r : , a ) - iQ„(r..,i - r.., a )].(27) 

With the help of this small source size (or large relative momentum) expansion, 
the two-particle Bose-Einstein correlation function can be written into the 
following Buda-Lund form: 

CAk, K) - 1 + X,(K) I'-W-q/, (28) 

|ff.(0)| 2 |G,(0)| 2 |7.(0;0)| 2 

This is the factorized Buda-Lund invariant decomposition of the BE correlation 
function; the resulting Fourier-transformed proper-time , space-time rapid- 
ity and transversal coordinate distributions can be of power-law, exponential, 
Gaussian or other types, corresponding to the underlying space-time structure 
of the particle emitting source: 

/>oo 

H*(Q=)= dr exp(iQ = r)iJ* (r), (29) 
Jo 

/oo 
dr7expHQ||T77)G*(77), (30) 
-OO 

h(Q:,Q..)= dTexp(-iQ.r : -iQ.r..)h(r : ,r..). (31) 



What are the fit parameters in eq. (|28|)? It is clear that apart from the shape 
parameters of the proper-time, space-time rapidity and the transverse distribu- 
tion of the production points two additional space-time parameters enter the 
fit: the mean proper-time of the particle production, r in eq. ( pQ) and the an- 
gular direction rj that enters the definition of Q= and Qn in eqs. (§])J(|). In turn, 
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the parameters (r, fj) are measurable from the detailed analysis of the multi- 
dimensional Bosc-Einstein correlation functions for any value of the rapidity 
and rat of the particle pair. The total longitudinal region corresponds to the 
region in (t,r z ) where particles of arbitrary momenta are emitted from. This 
region can be reconstructed in the Buda-Lund formalism, by determining (r , rj) 
and the widths of (r) and G* (77) at various fixed values of the momentum 
of the particles, reproducing the shaded region of Fig. 1 for each fixed value 
of the mean momentum of the pair. Such shaded regions are the same as local 
maps in cartography. If the momentum distribution of the produced particles 
is integrated over, the overlapping shaded regions are combined to a global 
picture of particle emission in space-time, similarly the way how local maps 
can be combined to an atlas in cartography by gradually displacing the centers 
of the local maps but keeping an overlap region between the neighbouring local 
pictures. 




Fig. 2. Space-time xdrjlure of particle emission as reconstructed from 
a combined analysis lilM of particle correlations and spectra in h + p 
reactions at CERN SPS in (t, r z ) plane. The momentum of the emitted 
particles is integrated over, from ref.Ej. 
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This programme has been carried out first by the NA22 collaboration in a 
Gaussian approximation, keeping the means and the variances only of the 
proper-time distribution, the space-time rapidity distribution and the trans- 
verse distribution coordinates of particle production. A locally thermalized, 
longitudinally expanding medium ( possibly corresponding to the vacuum filled 
by the sea of virtual quarks, anti-quarks and gluons) was found to be in a good 
agreement with the NA22 particle spectra and correlations, despite the fact 
that the mean multiplicity of the produced charged particles is only (n) = 8. 
See R. Hakobyan's contribution to this conference proceedings for the recon- 
structed space-time picture as well as for greater details. 

A similar reconstruction of the longitudinal space-time structure of particle 
production in 200 AGeV S + Pb reactions is reported in the contribution of 
Ster to this volume. 

For even further general possibilities about the structure of particle emis- 
sion, see refs. 



4 Earlier, Gaussian parameterizations of BE Correlations 

We briefly summarize here the Bersch-Pratt and the Yano-Koonin parameter- 
ization of the Bose - Einstein correlation functions, to point out some of their 
advantages as well as draw-backs and to form a basis for comparision. 



4-1 The Bertsch- Pratt parameterization 

The Bertsch-Pratt (BP) parameterization of Bose-Einstein correlation func- 
tions is one of the oldest, widely used parameterization called also as side-out- 
longitudinal decompositionEj'Efl of the correlation function. 

This directional decomposition was originally devised for extracting the 
contribution of the long duration of particles from a decaying Quark-Gluon 
Plasma, as happens in the mixture of a hadronic and a QGP phase if the 
rchadronization phase transition is a strong first order transition. 

The BP parameterization in a compact form reads as 

C(Ak, K) = 1 + A exp [-R 2 Q 2 - R 2 Q 2 - RfQf - 2R 2 ol Q l Q ] . (32) 

Here index o stands for out (and not the temporal direction), s for side and I 
for longitudinal. In a more detailed form, the mean momentum dependence of 
the various components are explicitly shown: 

C(Ak, K) = 1 + A(K) exp [- R 2 S (K)Q 2 S (K) - R 2 (K)Q 2 (K) 

-R 2 (K)Q 2 ~2R 2 0l (K)Q l Q (K)] , (33) 
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where the mean and the relative momenta are defined as 



K = 0.5(k!+k 2 ), (34) 

Ak = k!-k 2 , (35) 

Qi = k z .i - k z , 2 , (36) 

Q = Q (K) = AK-K/|K|, (37) 

Q S = Q S {K) = |AkxK|/|K|. (38) 

It is emphasized that the BP radius parameters are measuring lengths of ho- 
mogeneity, and in general characterize that region of space-time that emits 
particles with a given mean momentum K. Not only the radius parameters 
but also the decomposition of the relative momentum to the side and the out 
components depends on the mean momentum K. 

In arbitrary frame, Gaussian radius parameters can be defined, and some- 
times they are also referred to BP radii, when the spatial components of the 
relative momentum vector are taken as independent variables. It should be 
emphasized, that the BP radii that contain a well-defined mixture of the lon- 
gitudinal, temporal and transverse invariant radii. However, the BP radius 

ime where they 



parameters themselves are not invariant, they depend on the 
are evaluated □, reflecting space-time variances Era of the core 



of the particle 



emission. 



C c/,l (k, Ak) = 1 + A*(K) exp (-^.(KJAkjAk,-) , (39) 
A»(K) = [N C (K)/N(K)] 2 , (40) 
/.';.,.; K : = ((a* - - P s t)) - ((an - Pit)) B ((xj - 0jt)} e , (41) 

(f(x,k)) c = J dV(x,k)S c (x,k)/ J d 4 xS c (x,k), (42) 

where S c (x,Js) is the emission function that characterizes the central core. 
Note, that the tails of the emission function are dominated by the halo of 
long-lived resonancess Sh(x,\s.) and even a small admixture of e.g. r\ and 77' 
mesons increases drastically the space-time variances of particle production, 
and makes the interpretation of the BP radii in terms-jOf the total emission 
function S — S c + unreliable, as pointed out in ref. E3. 

It was shown already in ref . Ej that the duration of the particle emission 
contributes predominantly Aa the out direction, if the Longitudinal Center of 
Mass System (LCMS, ref. liil) is selected for the determination of the radius 
parameters. In LCMS, the mean momentum of the pair has no longitudinal 
component. In this frame, the BP radii have a particularly simple form, if the 
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coupling between the r x and the t coordinates is also negligible, (f x t) = (f x )(t): 

(3 2 {P) C (43) 

(44) 
(45) 

- /M)>c, (46) 

where x = x — (x). From the above, the advantage of the LCMS frame is clear: 
In this LCMS frame, information on the temporal scale couples only to the 
out direction and it enters both the out radius component and the out-long 
cross-term. Note that for cylindrically symmetric sources other possible cross- 
terms, e.g. the side-out or the side-long cross terms_were shown to vanish in a 
Gaussian approximation, due to symmetry reasons uu. 

An advantage of the BP parameterization is that there are no kinematic 
constraints between the side, out and long components of the relative momenta, 
hence the BP radii are not too difficult to determine experimentally. 



R 2 out (K) = (f 2 x ) c - 

R 2 lde (K) = (f 2 y ) c 
Rf ong (K) = (f 2 z ) c 



4-2 The Yano-Koonin-Podgoretskii parameterization 



A covariant parameterization has been worked. o.ut for non-expanding sources 
by Yano, Koonin and Podgoretskii (YKP) ElEa. This parameterization was 
recently applied to expanding sources by the Regensburg groupEZrEa, by allow- 
ing the YKP radius and velocity parameters be momentum dependent. This 
parameterization reads as 



C(Q, K) = 1 + exp 



-RlQl - Rj(Qf - Ql) - (Rl + i?f) (Q • U) 2 ] , (47) 



(Note that in YKP index o refers to the time-like components). When the 
momentum dependence of the YKP radii is explicitly shown, this reads as 



C(Q,K) = 1+exp 



-R 2 ± (K)Ql- 

tax) 



flf(K)(Q2 



Ql) 



Rj(K))(Q-U(K))' 



(48) 



where th^fit parameter U(K.) is interpreted@@ as a four-velocity of a fluid- 
element c3. This YKP parameterization is introduced to create a diagonal 
Gaussian expression in the "rest frame of the fluid-element" . 

This form has an advantage as compared to the BP parameterization, 
namely that the three extracted YKP radius parameters, R±, R\t and i?o are 
invariant, independent of the frame where the analysis is performed, while 
transforms as a four-vector. 
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The practical price one has to pay for this advantage is that the kinematic 
region in the Qq, Qi, Q± space, where the parameters can be fitted may become 
rather small. This follows from the unequalities Qf nv > and Qq > 0, which 
yields 

< Ql < Ql + Qi, (49) 

and the narrowing of the regions in Qq — Ql with decreasing Q± makes the 
experimental determination of the YKP parameters difficult, especially when 
the analysis is performed far away from the LCMS rapidities [or more precisely 
from the frame where C/ M = (1, 0, 0, 0) ] . Hence in practice the YKP parameters 
can be well determined in the LCMS frame, where the longitudinal component 
of the U is generally small. But in the LCMS, the interpretation of the BP 
radii is also simple, similarly to that of the YKP radii. 

Theoretical problems with the YKP parameterization are explained be- 
low, a) The YKP radii contain components proportional to that lead to 
divergent terms for particles with very—lnssi p t l3'C3. b) The YKP radii are not 
even defined for all Gaussian sources EZrE3. Especially, for opaque sources or 
for expaning shells with (f^) < (f*) the algebraic relations defining the YKP 
"velocity" para meter become ill-defined and result in imaginary values of the 
YKP "velocity" ,t3o. c) The YKP "velocity" is defined inlerms of space-time 
variances at fixed mean momentum of the particle pairs LxEij. Thus, for ex- 
panding systems, the proper interpretation of the parameter U is not a flow 
velocity of a source element, as thought before, but a combination of space- 
time variances of the source at a fixed mean momentum K . (Note, however, 
that for static, non-expanding sources the interpretation of UJL as the velocity 
of a Gaussian source can be preserved corresponding to refs.eJ'ES). 

In kinetic theory that provides the basis for hydrodynamics, the flow ve- 
locity can be locally defined as a weighted average of particle momenta, all 
particles being in the same cell in coordinate space. The local flow velocity 
u^lx) hence becomes a function of the position x but the momentum of the 
particles was averaged out, hence u¥-lx) is formally independent of the momen- 
tum. The four-current is defined cB as an average over the local momentum 
distribution 

/fl 3 k 
-yr^f(x,k) (50) 

where f(x, k) stands for the phase-space distribution function. The local flow 
velocity can be defined with this current as the unit vector proportional to the 
current: 

u"(x) = r(x)/^3^x)j^x). (51) 
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On the other hand, the YKP parameter U corresponds to a weighted av- 
erage of particle coordinates, all particles being characterized with the momen- 
tum K. Hence in general U(K) ^ u(x). The local velocity of particles at a fixed 
momentum is independent of the density distribution of particle production in 
coordinate-space. In fact, the only four-velocity that can be uniquely assigned 
to a set of particles each with momentum K is simply — {Ef^/m 1 "K./m), 
where m stands for the mass of the particles. Hence the YKP parameter U (K) 
should not be interpreted as a mean flow velocity of a fluid clement that emits 
particles with momentum K, at least not in the well-defined sense of kinetic 
theory. 



4-3 Comments on the hydro model parameterization 

The Buda-Lund hydro parameterization is recapitulated in the contributions 
of Ster and Hakobyan EirEj. The main idea behind the Buda-Lund parameteri- 
zation is to characterize the hydro fields (temperature profile, transversal flow, 
density distribution) with thcui means and variances only. Similar models are 
studied by the RegensburgLmEa as well as the Kiev groupsEJ. 

Note that n^(K) normal- vector rests in the (t,r z ) plane for the Buda- 
Lund parameterization. On the other hand, the C/ M (K) YKP velocity was 
thought to be some sort of flow velocity, that characterizes some local mo- 
mentum distribution. Hence, n M and U (K) are defined in different spaces: in 
the coordinate space and in the momentum space, respectively. In case of the 
Buda-Lund form, the coordinate-space interpretation of n is needed to obtain 
the expansion of the plane- wave exp(iAfc • Ax) in eq. (|2^), which is essential 
in expressing the correlation function in terms of Fourier-transformed proper- 
time distributions and space-time rapidity distributions in eq. (|28|). If the 
space-time interpretation of the Buda-Lund direction n(K) is lost, it becomes 
impossible to reconstruct the space-time picture of particle emission for sys- 
tems with strong longitudinal expansion. That could be the reason why such 
a reconstruction was not yet achieved by the NA49 experiment, that applied 
the YKP parameterization. However, the space-time picture of longitudinally 
expanding particle emitting sources was reconstructed from the detailed fitting 
to NA22 and NA44 particle correlations and spectra with the help of the Buda- 
Lund parameterization, properly preserving the interpetation of of n(K.) as 
a spatial direction. See Fig. 2, and thp_cjmtributions of Hakobyan, Ster and 
Seyboth to this conference proceedingsE3 ! M'E3. 
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5 Separation of particle sources in Buda-Lund type hydro models 

We discuss some results relating the separation of effective sources for identical 
and non-identical particles, even if both kind of particles appear from the same 
system that is assumed to be in local thermal equilibrium. 



5.1 Separation of effective sources for non-identical particles 

Recently, Lednicky and collaborators suggested to study non-identical particle- 
correlations, to lcaxn which particles are emitted earlier and which particles 
were emitted laterLj. The analysis of E877 data resulted in an effective sepa- 
ration of pion and proton sources at forward rapidities at the AGSLl 

Such separation of pion and proton source in space-time occurs as a natural 
result in the Buda-Lund hydro models and its various re-incarnations and 
improved modifications, because heavier particles are more frozen to the flow 
than the lighter ones. 

In the Buda-Lund parameterization, we have a K dependent normalvector, 
pointing towards the center of the particle production. 

In the LCMS, ref0 , we find 

n = n(x(K)) = (cosh[?j], 0, 0, sinh[r?]) (52) 

rj = rj(K) = V " ~ V , (53) 
1 + A?] z m t /lo 

where To is the central freeze-out temperature at the mean freeze-out time, 
and Ar] characterizes the finite longitudinal size of the expanding hydro source 
in space-time rapidities. 

For non-identical particle correlations^, e.g. tt and p, the velocities of the 
particles must be similar. 

If v v = v p , then we have << mf. For pions and protons with the same 
velocity, using = 140 MeV, m p = 1 GeV, T = 140 MeV and An 2 = 2 the 
Buda-Lund model yields: 



_ p — 7r y-yo /r,N 

V ~ V - — g — (54) 

where yo is the mid-rapidity, y is the rapidity of the pion and the proton, and 
3 is the numerically estimated coefficient. 

The spatial separation between protons and pions is about L ~ 10 fm/c 
at y — yo ~ 2E3. This can be used to estimate the mean freeze-out proper-time 
of particle production as r s = Lj {rf s — rjg) ~ 15 fm/c. 

1G 



5.2 Example for a non-Gaussian Buda-Lund correlation function: separation 
of the effective source for identical pions 

In the Buda-Lund type hydro models, the emission function S(x, k) is writ- 
ten as a product of the local phase-space density f(x,k) and a Cooper- Frye 
pre-factor d 4 T, fJ, (x)k fJi that yields the flux of particles through the_freeze-out 
hypersurface, or through a distribution of freeze-out hypersurfacesEJ: 

S(x,k)d 4 x = f(x,k)d 4 Y,-k (55) 



f(x,k) = 



(27T) 3 

exp 



k ■ u(x) — n{x) 



(56) 



T{x) 

d i T,(x) ■ k = m t cosh[r/ — y]H*(r) drfdr]dr x dr y (57) 

where g stands for the degeneracy factor, and s = —1, 0, +1 for Bosc, Boltz- 
mann or Fermi statistics, and an approximate boost-invariant shape of the 
freeze-out hypersurface distribution is assumed. Using the exponential form of 
the cosh[r/ — y] factor, the effective emission function S*(x,k) = =S(x,k) can 
be written as a sum of two components: 

S(x,k) = 0.5[S+(x,k) + S-{x,k)], (58) 
S+(x, k) = m t exp[+?7 - y]H*(r)f(x, k), (59) 
S-{x,k) =m t exp[-r 1 + y]H4T)f(x ) k). (60) 

These effective emission function components are subject to Fourier - trans- 
formation in the Buda-Lund approach. In an improved saddle-point approx- 
imation, the two components S+(x,k) and S—(x,k) can be Fourier - trans- 
formed independently, finding the separate maxima (saddle point) x+ and a;_ 
of S+(x,k) and S—(x,k), and repeating the saddle-point calculation for the 
two components separately. As a result, one gains the following non - Gaus- 
sian correlation function from the Buda-Lund hydro model specified in ref.EI 

C(Q=,Q h Q..,Q : ) = l + A*fi(Q||) exp(-Qj}iZf - QlRl - Q 2 ± R\), (61) 
fl(Qy) = [cos 2 (Q || A77) + sin 2 (Q || A?/) tanh 2 (77)] , (62) 

llt cosh[jj]. (63) 



Arj 2 An 2 T 

This result goes,beyond the single Gaussian version of the saddle-point calcu- 
lations of ref. LH3. This results goes also beyond the results obtainable in the 
YKP or the BP parameterizations. In principle, the improved saddle-point 
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calculation gives more accurate analytic results than the numerical evalua- 
tion of space-time variances, as it keeps more information on the shape of the 
correlation function. 

Although the above result is non-Gaussian, because the factor 
results in oscillations of the correlator, the result is still explicitely boost- 
invariant. 

Note that the oscillations are typically small and the Gaussian remains a 
good approximation to eq. (|6l|), but with modified radius parameters. 

The oscillations are due to a possible separation of the pion source to 
two components, due to an identical splitting of the Cooper-Frye pre-factor 
or the flux term in the emission function. We obtain two separate effective 
sources that create oscillations in the intensity correlation function, similarly 
to the oscillations in the intensity correlations of photons from binary stars in 
stellar astronomy. However, these oscillations in high energy physics are much 
smaller, as the effective separation between the particle sources for identical 
pions, rj + — rj_ smaller, than Arj, the width of the G*+(t/ + ) and the G*_(?7_) 
distributions. In stellar astronomy, the separation between the binary stars is 
much larger than the diameter of the stars, hence the intensity oscillations in 
the two-photon correlation function are stronger then in the Buda-Lund type 
hydro models. 

6 Highlights: 

The invariant Buda-Lund notation is introduced to describe Bose-Einstein cor- 
relation functions. This invariant notation scheme yields simple expressions 
not only for Gaussian but also for non-Gaussian expanding sources as well. It 
seems that the Buda-Lund form is the simplest and the most compact charac- 
terization of the two-particle Bose-Einstein correlation functions for relativis- 
tic, expaning systems. 

With the help of the Buda-Lund formulation, and a combined analysis 
of particle correlations and spectra Id , the space-time picture of the particle 
production in the longitudinal (t,r z ) plane, nan be reconstructed. Examples 
for such a reconstruction are shown in refs. rW™. 

We pointed out how to estimate the mean freeze-out time using non- 
identical particle correlations data, the Buda-Lund hydro model and the natu- 
ral separation of forward moving protons and pions in Buda-Lund type hydro 
models. 

Finally we have shown how Buda-Lund type hydro models can be rewritten 
to an effective, two-components sources, by splitting the flux terms. We found 
small damped oscillations in the intensity correlation function, reminiscent to 
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the intensity correlations of photons from binary stars. 
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